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ABSTRACT 


The use of low thrust devices with continuous and discontinu-~ 
ous thrust programs is investigated to determine whether or nota 
discontinuous thrust program will provide a greater payload and 
structure mass fraction, where the program must impart a specified 
amount of energy to the satellite in a specified total time, starting 
from a circular parking orbit. A discontinuous thrust program is 
developed, based on the elliptic orbit, using a perturbation analysis, 
and a series solution is obtained which permits investigation of this 
program up to angles of one radian either side of the perigee of the 
elliptic orbit. These results are compared with a continuous thrust 
program which gave a spiral orbit. Under certain conditions, where 
storage batteries must be carried as part of the payload but are 
available for use during the thrust program, the discontinuous thrust 
program is found to provide a greater payload and structure mass 
fraction. Further investigation at angles greater than one radian 


either side of perigee appears to be warranted. 
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LIST OF SYMBOLS 


gravity at surface of central mass 

radius at surface of central mass 
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angle of position vector in the plane of motion 
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non-dimensional time 

non-dimensional reciprocal radius 
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acceleration due to thrust 

acceleration ratio 
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mass of batteries 

specific mass of batteries 


specific mass of power supply 





ratio of we to M. 
exhaust velocity 

power supplied 

vector notation 

derivative with respect to 6 
derivative with respect to time 
perigee restriction 

Lagrange multiplier 

specific power 

initial cutoff angle 

number of thrust cycles 
factored ratio of “ie 

ratio of Mot to M 


ft g 
energy storage factor 





I INTRODUCTION AND SUMMARY 


Once in a parking orbit, low thrust devices may be used to 
alter the path of a satellite. In particular, it is anticipated that one 
use of the low thrust device will be to add energy to the satellite for 
an escape or a near escape condition. The exact path and thrust 
program must be determined by the mission requirement, but in any 
event, the system should be optimized for efficient launching. 

The problem of optimizing an arbitrary mission is a for- 
midable task to say the least. To simplify matters, it is assumed 
that the mission is to impart a specified energy in a specified total 
time to the maximum fraction of payload and structure mass. Todo 
this, one must make efficient use of the power supply and propellant. 

The usual variation procedures may be used to derive the 
differential equations for optimum direction and magnitude of thrust 
when continuous thrust is employed, as shown by Irving in reference 
(1). He has shown that a constant acceleration is optimum for the 
gravity-free case. No analytic solution has been found for the dif- 
ferential equations in the central force field, but Irving has made a 
numerical study and found the acceleration to be practically constant 
for various escape programs. 

In reference (2), Casey has shown the distinct advantage of 
energy addition by thrusting near the perigee of an elliptic orbit in 
such a manner that the perigee remains constant. His thrust pro- 
gram is formulated to require that the perigee distance remain con- 


stant throughout the thrust cycle, and results are obtained by 





ee 
numerical procedures. The reason for thrusting near the perigee is 
to take advantage of the greater velocity, since in any elliptic orbit 
the velocity is greatest at the perigee and since the energy addition 
rate is proportional to the velocity. Further, if the perigee distance 
is kept constant, then the perigee velocity will increase as the eccen-=- 
tricity of the elliptic orbit increases. 

As pointed out by Casey, there will be a period of coasting 
between each thrust cycle which will make a considerable contribu- 
tion to the total time. For this reason, we cannotsaysimply by in- 
spection whether the continuous or discontinuous method of energy 
addition will provide a greater fraction of payload and structure 
mass, where both methods are required to impart a specified energy 
in a specified total time. 

To investigate this problem, we shall use the idea of re- 
stricting the perigee distance, but only at the beginning and end of 
each thrust cycle. This permits added freedom in an effort to gain 
more energy. Since we are only considering low thrust devices, a 
perturbation technique is used with a simple variation to develop the 
first perturbation integral equations for the maximum change of en- 
ergy during one thrust cycle and for the Lagrange multiplier neces- 
sary to satisfy the above perigee restriction. The equations apply 
for any period of thrusting up to one complete circuit of 2w radians, 
but could not be integrated in closed form. By making a series ex-= 
pansion of the integrands, the equations may be integrated term by 


term to obtain results useful to an angle of one radian either side of 





aes 
the perigee. While this does not permit a complete comparison of 
these two methods of energy addition, it is possible to draw definite 
conclusions for angles less than one radian. 

After obtaining equations for one cycle of energy addition, 
the effect of a large number of cycles is obtained by integrating the 
changes of the various elliptic orbit parameters. 

To compare these results, an approximate solution for the 
spiral orbit, using alow thrust device, is obtained. Then a compari- 
son is made of the payload and structure mass fraction, where both 
methods are required to fulfill the mission of imparting a specified 
energy in a specified total time. 

Finally, since the discontinuous thrust program could 
generate energy while coasting, the use of storage batteries to aug- 
ment the power available during the thrust cycle will be considered. 
In one case, the batteries will be considered as part of the power 
supply; and in another, they will be taken as part of the payload, but 
available for use during the thrusting program. 

When thrust is only permitted up to an angle of one radian 
either side of perigee, the discontinuous thrust program was found 
to have no payload and structure mass fraction advantage unless bat- 
teries are employed. When the silver-zinc battery is considered as 
part of the power supply, there is an area of marginal advantage due 
to the low specific mass of this battery, but its poor recharging re- 
liability rules out its present use in such a system. 

The nickel-cadmium storage battery, which is very relia-= 


ble, provided an area of definite advantage for the discontinuous 
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thrust program when the battery is considered as part of the payload. 
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We shall first write the equations of motion for the satellite, 
making the usual simplifying assumptions of a very large, central, 
spherical mass and take the coordinate system at the center. We 
shall further assume motion ina plane, keeping the thrust in the 
plane of the orbit, and neglect the effects of other bodies. Then, 
using plane polar coordinates, as in fig. 1, and the dot notation for 


time derivative, we have: 


Fim = OT (2.4) 
tm 2 TR +jne —- 
- ee: é Fr A 

fake = ( (J ioe ag ice es ef | (2255) 


Equating equation 2.3 to the acceleration due to thrust, Q, and grav- 
ity, we have: P 
— 2. oy (7 A Nn? 8} = C(a cos( B+) oe Re ) 
i(n-né*)+ I (Fat ne 
+ j (a sin( B+) ) (2. 4) 


where R, is the radius of the central mass, and 9, is the acceler- 
ation of gravity at n=R,. Since 9/4e is small for a low thrust 
device, we shall first take Q=0 and obtain the usual solution in 
dimensionless form. The thrust will be considered later as a per- 
turbation. Choosing © as the independent variable, and letting ( )! 


be differentiation with respect to © , then with U as the dependent 
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variable and h as the dimensionless angular momentum, we have: 


R, 


U A (2. 5) 


ne = a/R. Roh (2. 6) 


but when Q=O , 4 (ne )=0O ; therefore hi = constant = 


Woe 
C= V9/Re hu" — 


A =~Ha.R, hou (2. 8) 


oe ce 2 “0 
IT al 90h U U (2.9) 


12 2 3 
ae, om Go Ay U (2. 10) 


Substituting into equation 2.4, taking @Q=0O » and simplifying, 
we have: 
ty | 


UW + U = a (220) 


9 


A solution may be written 
S 
oe Ue = yall + @, cos (@ - @p)| (2, 12) 


where €, is the eccentricity and ©p is the angle of the perigee 
position. The non-dimensional time, AT , for any portion of the 


orbit is given by: 
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ey 0, 6, ; 
= 4 = fo = 4e de = 3 6 
q JEG u je | a HR) 6 h (1+ e,c05 @)* (Zain 
@, - 4 


using equations 2.7 and 2.12. The integral can be evaluated to give: 


= 
: [ +, 


O 
2 [e SE Pa fal Coll-C2 tan 2 
Beye) Pie 148) — 7 esses 
(1-e2) a (ite) +Cl €,) taW 2 


In equations 2. 14 through 2, 26, we shall essentially follow the 
work of Irving, reference (1), since his approach provides a conveni- 
ent parameter for computing the mass fraction of the satellite. 

The mass of the power supply, Mw , is assumed to be pro- 


portional to the power supplied, P 
Mw = KP (2. 14) 


where | _ is the specific mass of the power supply. Itis further 
assumed that the maximum power is utilized when thrusting, appear- 


ing as kinetic energy in the exhaust: 
= : 4 
esl G oC (2. 15) 
where C is the exhaust velocity, ”™ the mass of the vehicle at 


time t , and -™ is the propellant flow rate. We also have the 


acceleration of the satellite due to reaction: 








a = thrust a mec (eo) 
mass wom 


Using equations 2.15 and 2.16, we eliminate C 
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2P ~ wy 0? i 5 ae (Zak) 


Equation 2.17 may be integrated to obtain an expression for mass as 


a function of time: 





| Taeee oe At 
O 
where Mp, is the initial mass of the satellite. Ifwe let Mp be 
the propellant mass, then the ''burn out''mass, M, , at time 7, 


iS given by:. 


Mi — M. _ M, (2, 19) 


Substituting in equation 2.18 and using 2.14, we have: 


Me) 4, K Me i (2. 20) 
pal INA 


To optimize the powersupply, we write the relation between 
the payload mass, M, , the structure mass, Mes » and the cther 


masses. 
aM RS ocala al (Zea) 


Further, let 
ie 
ye - > \ Canie (eae) 


Then, manipulating equations 2.19, 2.20, 2.21, and 2.22, we have: 


MitMs Mw fot 1 2, 23 
Mie / nae CM Meyer Ve 





aon 
There is an infinite number of programs for Q which will keep Y¥° 
constant while fulfilling the mission as ie jhe is varied. Thus, 


taking the derivative of equation 2. 23 equal to zero, we have the 


maximum of (M,+M,)/ Me with respect to My/M. when: 
M 
iG =y -y2 (az 


Substituting equation 2.24 into 2. 23, we have: 


Nees Geayeye (225) 


o 


Using equations 2,21, 2.24, and 2.25, we have: 
eee, — 8 (2526) 


From equation 2.26 we see the range of YY must be zero to one, 
since Mp/Mo Must 16 in this range, Further, tor a paviiculan 
mission, we see from equation 2. 25 that the structure and payload 
mass fraction will be greatest when Y is aminimum. After de- 
termining the minimum value for Y¥ , the optimum power supply 
mass 1s given by equation 2, 24, 

We are now ready to investigate equation 2.22. In reference 
(1), Irving has shown analytically that the acceleration must vary 
linearly with time in gravity-free space, if Y is to be a minimum. 
Further, when the mission of specified energy in a specified time is 
considered, the optimum acceleration becomes constant. For the 
case of a central force field, Irving made a numerical study using 
continuous thrust, and found the optimum acceleration to be practi- 


cally constant. Tnus, the spiral orbit with constant acceleration is a 





<i Oe 
"near optimum" solution which will provide a fair comparison be- 
tween continuous and discontinuous thrust programming. 

If we consider a discontinuous thrust program applied to an 
elliptic orbit, it seems clear that the thrust should be applied near 
the perigee, where the velocity is greatest, to impart the greatest 
energy. If the thrust is parallel to the velocity vector, we will im- 
part the maximum energy during any given thrust cycle. Such a pro- 
gram will allow the perigee distance to increase, and with this in- 
Geease there is a decrease of perigee velocity, thus decreasing the 
energy imparted during a subsequent given thrust period. To obtain 
the advantage of thrusting near the perigee, but still keeping the 
same distance for a velocity advantage in subsequent cycles, a pro- 
gram must be found that imparts maximum energy during any given 
thrust cycle while providing no net change in the perigee distance. 
After obtaining the restricted perigee equations for one cycle, the 
programming of thrust time will be considered and’ Y will be com- 
puted for a large number of cycles. 

After developing the restricted perigee equations, we shall 
obtain an approximate solution for the spiral orbit with continuous 
thrust to compare with the restricted perigee thrust program. Fi- 
nally, we shall consider the use of batteries with the restricted peri- 
gee thrust program. The batteries can be charged while coasting to 
provide added power when thrusting. In one case, the batteries will 
be considered as part of the power supply, and in another case they 


will be considered as part of the payload. 


eS 
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When Q@ is not zero, we have a small non-dimensional ac- 


celeration: 
2. 
os 


Go 


since we are considering only low thrust devices. This suggests the 


<< ] (3, 1) 


use of a perturbation analysis for each cycle. 


Since # is no longer constant, we must rewrite equation 


29: 
A ee Got ye eee) (3, 2) 


anaeor the i component of equation 2. 3 we also have: 


Ld (nt) = gah 2.3) 


By using fig. 1 and simple trigonometry, we write: 


—~Uanp + Ucosp 





SIN(P+B) = yaerue (3. 4) 
and 
Lr‘eos b - U sis 3. 5) 


cos(p+B)= fara 
Sustiiuting equations 2.5, 2.10, 3.1, 3.2, 3.3, 3.4, and 3.5 into 


equation 2.4, we have the radial equation: 


Nee f 


“ fu | Pe ‘cos d + USING 
men) ve eat tour) 


and the tangential equation 


(3, 6) 








ae 


gag -uU’sind + Ucosdh 
- oe Tur +u? 


[Sia) 


/ 
We get some simplification in equation 3.6 by substituting for fh 


Jo ae Vute wu’ sin 


ee Wee 


irom equation 3.7: 


For the perturbation solution, we take te constant and U, as 


given by equation 2, 12 to write the assumed solutions: 


U = Us + EU, + E*U2 + CPU, + --- (3.9) 
and 

fh= th, + Et, + C7, + O7h, Fo (3. 10) 
where U,=U;,=°-° = Y= ee oO) at ©, , the angle where 


the thrust for a particular cycle is started. Substituting equations 
3.9 and 3.10 into equations 3.7 and 3.8, we obtain the first perturba~ 


tion differential equations: 


-U, sin d + Us, cosh 


/ 
h, =< 3 2 Pe. (3. 11) 
tu; Wu: + US 


“ Ve a 5 
UV a hy 4 Vue tue yg (Sa02) 
h, he, U> 
For a particular cycle, where the thrust is started at angle ©, and 
9-08 =—27T , the solutions may be written: 
o 7 
h, (8) = f ne Cmicate (3. 13) 
Oo, 





alse 
and 


e Y ae 
U(@)= | si(6- aes \k (x)dx + Vue + UT et sp} AY 
eb 


o i 


(3, 14) 


sein 


For aunit mass, the energy is given by 


F = kinetic energy + potential energy 
" d 
a 2 calihal 
a g(RP + Reet) + GREY St 


rE @E (Ru 4 Bu?) oh F(t) 


ee: 
= goRe | | a oa | (3, 15) 


by using equations 2.5, 2.6, 2.7, 2.8, and 2,12, and where the po- 





tential energy is considered zero at n=R,. By substituting equa- 
tions 3.9 and 3.10 into equation 3.15, we obtain the first perturba- 
tion energy change: 


Aig 


=F (uetus)h, +(pu-iu, + Pus us (3. 16) 
E 4o = 4.Re 


The perigee distance ratio is obtained by evaluating equation 


Zeke at <) ae 


R ina 
= = J (Semler) 
Me Ae 


We could require that this quantity remain constant during the thrust 
cycle, as done by Casey in reference (2), but we should have a gain 
in the energy imparted if we require that the perigee distance ratio 
be equal berore and aiter the eycle. We shall pursue this latver 
method of perigee restriction and will demonstrate a slight gain of 


energy over the method of continuous perigee restriction. 





ine €ccentricity may be evaluated betore andeatier taemtanuse 


C= Jug? + Uh" (3. 18) 


by using equation 2,12. Then, substituting equations 3.9, 3.10, and 


ey cle: 


3,18 into equation 3.17 and equating the results before and after the 


thrust cycle, we obtain the first perturbation perigee restriction: 


Pe - zp ye. Uo So(I+ee) 
PR=O= 2[ hus +hue- 5 get | + 


foe, UES Wi <2 ta, WT 
(3. 19) 


evaluated at ©, , where ©, is the angle at the end of the thrust 


cycle. 





Sac 
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PERIGEE CYCLE 


To find the thrust angle which will impart maximum energy 
Curing the thrust cycle and still satisfy the perigee restriction given 
by equation 3.19, we shall take a simple variation, using the La-~ 
grange multiplier, equal to zero. We write the non-dimensional en- 
ergy, equation 3.16, plus the perigee restriction constraint, equation 
Slo, with the Lagrange multiplier, d , where A is a constant to 
be determined, and, assuming € = constant during the cycle, take a 


variation on p equal to zero: 
Ae 
———— + PR | =) 
S E ak, YA (4. 1) 


For convenient notation, we rewrite equations 3.16 and 3.19 as 





as An, te GU, + Cu; (4, 2) 
and 
PR= 0 = AUX, + Bu, + Cus (4, 3) 


/ 


where A, A , 6, and C€ are shown in equations 3.16 and 3.19. 


Here, we note,by taking the difference of equations 4. 2 and 4.3, that 


DE 
E Fo Ko 


is the energy change for the thrust cycle. Substituting equations 4. 2 


= (A-A’) H,(4&) (4, 4) 


and 4,3 into 4.1, we have: 
(A+ AA)SA, FU +L BSU, + CSu/ ] = 0 (4. 5) 
Before taking the variation, we must first integrate equation 3.14 by 


parts to obtain the more useful form: 





Sees 


o 
U,(0)-\ sie) ee ae Me sind + 2; eo pT RAY (4, 6) 


6, Qo 


Then, taking a derivative of equation 4.6 with respectto 96: 


8 
U6) Bowe: yO sia GSP - $3 | Ssin(6- -v) ty (pL) d¥ (4, 7) 
8, aa 6, 


Misme equation 3.11, we note that: 





Oz Oz 
/ STC Os - U,gs/nv 
$4,(6,)=5 SB/rdy = | ; . p {od Av (4. 8) 
Se 6 le yu; gS 


Now, taking the variation of equations 4.6 and 4.7, using 4.8, we 


substitute into — 4. 5 to obtain: 


seq. 








Us tls 


me + A) Vie 
— Z. tes ¥) +Ceos(0,- Wi) cosd - 





(a2 -|{- or sine at SING +U, ‘eo Soa ¥ 


(459) 


For equation 4.9 to be zero, we must require that the integrand be 


zero, so that 





al fe 


taf, 28 : : 
U, ain p) AEA + Se (cor(e.-¥)-I) m7 EE ING, -¥)] = 


es 





=— U, ‘cospf AtME + SF (c05(6,-¥)-1) = a SF) = 


ees 
oe We PMe (B sin lO, 1D) 2 aa Ga v)| (4, 10) 


Substituting the expressions for A A Ne P (s a ane C from 
equations 3.16, 3.19, 4.2, and 4.3, using equation 2.12, and taking 


Op =O , equation 4.10 becomes 


é 
SIN © Ae, (lt €&o) 
ia ? oe nes ee ar EVN 


where A is a constant to be determined. It will be shown that the 
shift of Op is negligible for the symmetrical thrust cycle, so that 
a choice of zero for Op presents no problem in going from one 
thrust cycle to the next. We now have the thrust direction program 
given by equation 4, 1l. 

Substituting the expressions for A and A’ and equations 


Z,12, 3.11, and 4.11 into 4.4, we have: 


AE 

ef ea.Rbe 

-{p (14e,) (14 e, cos [1+ OF +2 6, Cos¥ +2 AC (cos ¥- ] - €o(it@o) Neo SIM ey AY 
aS 


(1+, cos ¥)° Jit e2 +2C,CosP [1(1te Cos) “[i+e24 2e,CosV +2 AE, (cos Y- =p) Jo + 


tM ES (Ite) acnek ¥ 
(4, 12) 


Also, substituting equations 3.11, 4.6, 4.7, and 4.11 into 3.19, and 





-18— 
using 2.12, we have: 


PR=0= 


© 
fi PeLe, Fo? Y £2 (142, cosW Al 1 +05 +20, Cos ¥ #222, (cos -1}] (cos P- i) AY 
a 
Cit e,cosy)? Jite? +20, CosY yor aCOn'e) [lt@e+2eo cst Fares (cos ¥- Mie Se 
+ Ver (ite.y sin y 


(4. 13) 


From equation 4.13 we determine the value of the constant A@€@. 
for any particular cycle, since A always occurs with €, in the 
two equations 4. 12 and 4. 13. 

Since it is not possible to integrate the above equations exact-= 
ie we must use numerical methods to carry the investigation to 
large angles, This would be fairly simple for one cycle, but when 
one considers a large number of cycles, a computer program appears 
to be mandatory. We shall not pursue this route, but instead, we 


shall present an approximate evaluation which is useful to angles of 


one radian. 
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Vo APPROXIMATE EVALUATION BPOR-ONE CYCLE 


To make an approximate evaluation of the optimum restricted 
perigee equations for one cycle, we shall make series expansions and 


Mavegrate term by term. First, using equation 4.13, we may rewrite 





4,12 as: P 
AE pR (Huse se eth I+ e,cosw)@[ +e, +2 en Cos¥ +226, (cos ¥-/) | dV 
eo ee 
E4,Rohe ; DENOMINATOR 


(5. 1) 


since PR=O by the appropriate choice of  A@, and where 
DENOMINATOR is the denominator of the integrand in both equations 
4,12 and 4, 13. 


Now consider the function 
Ce. 

OF \ d a (5. 2) 
O, 


where 


: (i+ @,cosy)*[ (I tes+2€,cos¥)b +22, (cos¥-i)] +X eS (1 +e.) Siw? Y (5, 3) 
(1+@,cosy)? Fite? +26,cos¥ 


We notice that 





(5, 4) 


has an integrand identical with equation 4. 13, and 
0, O2 
JG} 2. J _ ® 
ab =3,\¢ | = | fay oa 


bel 6, b=] 0, as 








= 20s 


has an integrand identical with the non-vanishing term of equation 
5-1, since PR=O . Thus, if we perform the integration of equa~ 
tion 5.2, we mayevaluate A@. and the increase in energy by 
taking the derivatives indicated in equations 5, 4 and 5. Bs 

Expanding in series, we collect the coefficients of the various 
terms of d and integrate. Noting that the result is symmetrical 


about @G=O , we write 








a 
all b N , ae Ne ee N3 0, - ‘f 
ZG “a $a) | * (eb (te)tiz0 ( #2.)¢ 16, 080” (5. 6) 


where 


N,= 2@, +€,-27+C? (5. 7) 


Nine 160, +1905 -eF -e, -22 705 -327e,+27+4+20 772 + 


+3407, -47%41207-3 04 (5. 8) 


4. a 
N;= be Po Cis e* 4242 CL eer= 18 Be +2Co -34t9G TC. -2XRGWE TE, 
2 
~ 1746 Tec +344 7e, -¢0 4/272 refs zezer%eg +2802 Tees - 


S = Rew Fo 2 _ 
-922 070, + 3202 +1140 Fes +2640 772, -3007T4-4207 %e;, 


G 
-930T%e, + 8407%- 54075 + 907 ae 


and for convenience we have used 


ee 


et ess (5, 10) 





Now, taking the derivative of equation 5.6 as indicated in equation 5, 4 


equal to zero, we satisfy the perigee restriction of equation 3,19 and 





Spa 
solve for A@e,: 


Pees 42 30ise,'- 30726 -)065¢e°-9936e, -3497 au 
NE |S a 


ce Deen 
a (i+eé.)7 16,800 





Hiere, we note that equation 5.11 would be much more complicated if 
a non-symmetrical thrust program were used. However, since we 
can see from the form of equation 5.6 that this will have no energy 
advantage, a symmetrical thrust program is chosen to simplify the 
results as much as possible. 

Similarly, taking the derivative of equation 5.6 as indicated in 
equation 5.5, and using equation 5.11, we have the energy for one 


ait Of the cycle: 
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To determine the advantage of this method of perigee restriction 
over that used by Casey, the energy equation given in reference (2) 
was expanded in series and integrated to the same number of terms 
as we have in equation 5.12. The two energies are identical to the 
last term, but if we compare the last term of the constant perigee 


energy expression: 


~ 2 42 
F 2720€° +5190€2 +1250, - 2000 2 - 255005 - 42300, -2335 ra 
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with the last term of equation 5.12, we see that for values of €, 





=22= 
between zero and one, the constant perigee restriction yields less 
Srecney. 
To find the changes of the elliptic orbit parameters U, +h, 


and © , we note from equation 4. 4 that 


| _ aE 
(A-A’) Ego Re —_ 
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where the energy is given by equation 5.12. Substituting the expres- 
sions for A and A* and neglecting the higher perturbations since 


€ is small, we have 


4 AE 
mm = o Senha Bele 
Ah = E%,(0,) = — ; 9 Re ( ) 


Now, from the perigee restriction, we may write equation 3. 17 before 


and after the thrust cycle: 


[+ @, (+e, tae 


fe (e+ ahh st Sa 





imeem which we solve for the change in eccentricity: 





be AE 
se = 6 Gyan 70 Re (5.17) 


using equation 5.15. To see if there is a shift of perigee, we assume 
that the equation for U , at the end of the thrust cycle, has the form 


Pie 2: 
[1 +(€,+4€) Cos (62-4 Gp) 


ae 
Zz belo 
(B,+ak) =" 
where we have provided for a shift of the perigee with the term A Op ° 


imen, using equations. 2.12, 3.9, and 5.16, we can write 
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re } r(Gt2&)cos(6,- 6p) ~ (142, c080,) (1 + aon 
vel Go A (5.19) 
and 
L. 
a E{ G. tabd® (5. 20) 


Substituting these two equations into 4.3, with the expressions for 


A’, B , and C , we obtain the result 
CoS (O,- 462) -C0s@,) Cos 6, (51M (61 -48p)-Siv 6,) 5/18, =O + OE) (5, 21) 
2 P 2 


Thus Ap is at most of order € . Using this result, we may 
evaluate the terms of order € in equation 5.21, using 5,15 and 
5.17, and we find that they are zero, Thus, there is no shift of the 
perigee within the magnitude of the first perturbation as would be ex- 
pected with a per cetiean thrust programa, 

Before proceeding with the discontinuous thrust program, we 
shall obtain the approximate solution for the spiral orbit with thrust 
mimeeme LO the velocity vector. or this case, p is zero, and equa- 


tion 3. S$ becomes 


2 | 
i — eeu (55 22) 


From cquation 3. 7 we have 
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/ 
/ & uU\ye7T7B 
aae= [i+ (Z) | (55-23) 
Taking twice of equation 5. 23 equal to the derivative of 5. 22 and ex- 


panding, we obtain 


— Cl 
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uul=-zefie2 AL (4). Tu (5. 24) 


where we have assumed starting from a circular orbit with ( U/L) 
and (U‘/U) much less than one. Integratine fromm Ce couc amc 


have: 
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2 ies uUu 
Ut= uz- 466 -2((UU sees 
(9) 
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To determine approximately the small quantities not yet integrated, 
let us take for the first approximate solution 
(ae GC 
eee, = 4 20 (5.26) 
‘Wogemis merely the first result from equation 5.25. Then, taking 


derivatives of equation 5. 26, we have approximately: 
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Subetituine equations 5,27, 5:28, and 5. 29 into 5. 25; we obtain the 


result: 





Wig Ue +226 fe 7 a (5. 30) 
Ua U, 
where Ug, is given by equation 5.26. We may now use the second 


term of equation 5.30 to choose values of € and UU, , assuming 


Ul, is given, to insure that our results lie within a certain accept- 
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able error, We should also note that (€/U,) must be much less 
meat, One it equation 5. 30 is to have meaning, 

io determine the time, we use part ofvequatron 2 omimncecun= 
stitute equations 5, 22 and 5.30. Then, using equations 5. 26 and 5. 28, 


we obtain 
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Here, we see that for specific values of Ug and € , we will have 
eplanger error im e¢quation 5,50 than in 5.51, since {2° 1s lessitaan 
U,. and U, is less than one. 

Te determine the enérgy input, we subsiiture equations 5,42) 
5.27, 5.28, and 5,30 into 3.15, and take the difference between the 


imal and initial values: 


or; 7 2 Ue Ua) — e*( ce - tt.) | (5, 32) 








where again the error may be estimated by the second term. Of 
these three equations, 5.30, 5.31, and 5.32, equation 5.30 will have 


the largest error as Ua, becomes small, and it should be used for 
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choosing values of € and U, . 

We have reached the point where we must decide how to pro- 
gram the cycles for the restricted perigee. To do this, we shall 
choose a ratio of (Ro/rt) =0.95 tobe the maximum value for U . 
For the case of the earth, this corresponds to an altitude of 208 miles 
and is adequate to insure negligible atmospheric effects. 


If we form the quotient using equations 2.13 and 5.12, 
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specific power == --~ —= Sy oye. 
pecific pow V ER, AT (5. 33) 
then expanding AC in series and performing the division 
mie | © 2 h= (Seamer 
— a. a of > Pa oo a eee? (5, 34) 
oe L+@, G CS ae SEC 


is equation is plotted in fig. 2 for various values of “Ge. ) with 


corresponding values of (e / eo Re ) for ie = 0.95  , where 


BS is defined as 
(Peres 
a= pe Guest 
oO 


Also plotted in fig. 2 is the specific power for the spiral orbit for the 
same values of (E / Go Re) » We can see that wnen the eccen- 
tricity is small, there is not much difference between the two methods 
near the perigee, but the discontinuous thrust program becomes inef-~ 
ficient if we proceed too far from the perigee. This is because the 
thrust must be directed away from the velocity vector to satisfy the 
perigee restriction. Thus, for a given energy increase per cycle 
when the eccentricity is small, € must be larger with discontinuous 
thrust, since the thrust is not on for the full 277 radians. There- 


fore, 4Y is larger, and the payload and structure to initial mass 





se 
ratio is reduced, as seen in equation 2,25. If we manipulate equa- 


tions 2,13 and 3.15, we can obtain the time for 277 radians of trav- 


ein terms of the total energy; 


s 
_ ©yR'f,. — |] 2 
ae ele! ae (5. 36) 


Which holds for both elliptic and circular orbits. Now, if the mission 





is to be accomplished, both methods will have about the same number 
@meycles, where in the case of the continuous thrust method, cycle 
means C7) radians of travel. Thus, we want to impart about the 
same amount of energy during each corresponding cycle. 

We can also see in fic. 2 that tor larger yalucc @itccecn rier 
say 0.75, the specific power is much greater for the elliptic orbit. 
In this case, for agivenenergy increase per cycle, AY for the 


Seine! Orbit will be greater. 
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Vi, APPLICATION TO MULTI-IEE Oo veren. 


To find the total energy input, total time, and Y , we must 
make a summation of the results for each cycle. However, since we 
are assuming a large number of cycles, with small changes for each 
cycle, we may write integrals for the summation operation and ob- 
tain our results by integration. The thrusting time, Mt , is taken 
as constant ior all of the cycles, and the variable of integration ts 
de 

We must first express the cutoti angle, (OG gyim rermianomee 
because it must be varied slightly to keep the thrusting time constant. 


Expanding equation 2,13 in series, we have, for the symmetrical 
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fo 6 Oe i 6 56 CCnee 7 ae (6, 1) 
‘he Pico] 3 0+e) *  Gocieer® 
Then, if we let 
ie /2 
‘fe eee 
R R, Zz (6. 2) 





we can write approximately 


= Ayre fi-ge@ eS 8 (6.3) 


where k is the initial value of Oo, - Care must be used in 
choosing R sothat ©, islessthanoneat C=] j; otherwise, 


the series expansions will not be good approximations. 
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Substituting equation 6,3 into 5.17, using 5.12, and replacing 


the finite ''delta'' quantities with differentials, we obtain 





. PE ge) 12.€°-29e°-32e-Il 44 
. 4.% ( ie 360 Ce nke Ryde (6.4) 


The term di is in effect one, that is, one cycle, and (de/E) is 
tae srnall change of eccentricity during that one cycle. To periorm 


the integration, we shall take € as constant to obtain: 
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fi py 3821+1828E-577e"+216e7 * 
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where ©; is the initial eccentricity and @¢ is the final. N is 
the number of thrust cycles, and the total thrusting time is given by: 
ee =e eae (6. 6) 
We find the energy input very simply by replacing the finite ''deltal 
quantities in equation 5.17 with differentials, using 5.35, and inte- 


grating to obtain 
(Ee, - €;) (6.7) 
Then, using equations 2.22, 3.1, and 6.6, 


yea Sg, Nat = K gerJEnat (6. 8) 


To find the total travel time for C¢7YF radians, we write equation 
5.36 in dimensionless form, replacing the finite "delta' quantities 


with differentials; then substituting 3.15 and 5.35, we have 
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dU = (6. 9) 
where again, adi is in effect one. Substituting equation 6. 4 for 


di and neglecting "end effects" in the integration, we have 


“ o 
= 12- 40e-72e% 
ue +X paswie Sh (6. 10) 
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The "end effects" we mention are the initial cutoff angle and the final 
cutoff angle not being equal so that we have a small amount of time 
merelected for each cycle due to the shift of the cutom angle; the dima 
energy is added before starting the final coasting orbit; and the initial 
energy is added during the latter part of the initial coasting orbit. 
This might be clearer if we note that the integration gives us 


ae oe) ee ee ee Goi 
is ei Se, Ce, es ( ) 


where A Te, is the travel time for 27f radians at eccentricity 


Ci > whereas the true mission time would be given by: 


—_— —_~ —_ —~ — ere were iia if +a Tl (6, 12) 
C=AT +aT tdt ral, + d,,t a0 ees a 
tft E*2 
where aT; is the thrust time for the initial cycle, 
C@-=e cay aoe , and o iF is the cutoii angle shitt 
a Cr J 


time. If the number of cycles is not large, then equation 6. 10 should 


be corrected for ''end effects'’ by adding 


Cx ==) AT. - AT é 
/ = = - 13 
CORR ck (| oe c. As (Ons) 
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where @ is the initial cutoff angle given by equation 6, 2, ye is 
obtained by evaluating equation 6.3 at @=@r . The cutoff angle 
shift time correction is an average value obtained by using the first 


term of equation 6,1. The other corrections follow simply by com-~ 


paring equations 6,11 and 6.12. 
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VI. RESTRICTED PERIGEE CYCLE WITH 


BATTERY STORAGE 


We must still consider the possibility of different values of 
specific power plant mass, K , for the two systems, since we have 
the capability of storing energy with the discontinuous method: while 
coasting if we provide batteries, The energy generated by the power 
supply while coasting could be stored in batteries and then used in 
addition to the power supply while thrusting. 

To study this possibility, we shall first assume that a battery 
pack is included with the power supply and that other conditions of 
the thrust program are unchanged; later the battery will be consid- 
ered as part of the payload. We shall further assume that during the 
early part of the program, when € is small, that the coasting time 
is not sufficient to fully charge the batteries, butas @ increases, 
a point is reached where the batteries do become fully charged, and 
are fully recharged during each subsequent coasting period. We 
have equation 2.18, evaluated at the "burnout condition": 


Tb oy 2 
pies fy ae ~{ Mw Oo dt 
0 


eM 2 KP (7. 1) 


where we have multiplied and divided by iS icy) , and must now 
consider FP asa function of time. When the batteries are not fully 


charged, the power is given by 


ne Ma Aan Wave 
a. IK Pete a2) 


neglecting battery losses while charging, where Mg is the mass of 
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Vit, RESTRICTED PERIGEE @ 2 be ayia 


BATTER Gs ler. cE 


We must still consider the possibility of different values of 
specific power plant mass, K , for the two systems, since we have 
the capability of storing energy with the discontinuous method: while 
coasting if we provide batteries. The energy generated by the power 
supply while coasting could be stored in batteries and then used in 
addition to the power supply while thrusting. 

To study this possibility, we shall first assume that a battery 
pack is aanaed with the power supply and that other conditions of 
the thrust program are unchanged; later the battery will be consid- 
ered as part of the payload. We shall further assume that during the 
early part of the program, when @€ is small, that the coasting time 
is not sufficient to fully charge the batteries, but as @€@ increases, 
a point is reached where the batteries do become fully charged, and 
are fully recharged during each subsequent coasting period. We 
have equation 2.18, evaluated at the "burnout condition": 

Mo _ » Me Khe OF yy (7. 1) 
M, ie a a P ; 
where we have multiplied and divided by GOIinoy, , and must now 
consider F asa function of time. When the batteries are not fully 
charged, the power is given by 
P= Ma A Teyere 
IK ee 2) 


neglecting battery losses while charging, where M2 is the mass of 
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the electrical generating unit, K is the specific mass of the gener- 


Semne unit, A 


Te 


is the constant thrusting time per cycle, and 


eyele is the travel time for ¢W radians given by equation 
5.36. The batteries will be fully charged when 
M 
M batt. _ g (aT, ™<¢s av) ms, 
kK kK Orbit 
batt. 
Or 
AT ..2 are M pate. Ko an 
orbit M IX 
d batt. 
where il ee and |x 


batt, 27e the mass and the specific mass of 


the battery, respectively. From equation 7.1, we define me 


, cor- 
responding to 2, 22: 
T 5 7 
ye. Bl Me Ope Ee ae 
2 K P es p (7. 5) 
o o 
leap since © is constant, and My. = Me + Mapai ° 
eg? Ms (+ Fgh) FO) 
‘acs = (Tate) 


where 7 os) is introduced to remind us that the integrand is 


zero when not thrusting. Using equation 7.4in 7.2, we can write 


the maximum power available when the batteries are fully charged: 


M att. IK 
— Me 22 Je pone eee ee 
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Bein? €quations 7.2iand (2.7 a fom 
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CRIT 
where t rit is the time at which the batteries become fully 


charged. Now, the integrals in equation 7.8 may be treated as sum- 
mations since the thrust is on for the constant time AT during 


each cycle: 


Nergir — 
ergi Mrat\ VR. > (N=Newir) AT 
yes K (| + pai) 3: (AT) TE *( jp sluiiaieaten (7. 9) 
d ee a AT M4. K batt. 





where VJ is the cycle mumper aiter Whientine Dareriecrere 


ferme eo 


fully charged. Using equations 6.9 and 7.4, we can solve for the ec- 


centricity when the batteries become fully charged: 


| 277 
S| -5| ____ (7. 10) 
: SL AT + Det 
M4 Kyat. 


Then, using equations 6.4 and 6.9: 
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Integrating this equation and dete rmining N-N from equation 


Crit. 
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6.5, using 7.10, equation 7.9 ae 
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Now if we form the ratio of equation 2. 22 for the spiral orbit and 


equation 7.12, multiplied by the factor J} +H » we have 
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When F Rye is less than one, the spiral method of energy addi- 
tion must have a Y less than the Y for the restricted perigee 
method, and by equation 2. 25 the spiral method must have a greater 
mevioad. If FR y2 is greater than one, then there may pewareoce 
sibility of carrying a greater payload by using the restricted perigee 
method. This will depend on the ratio of K/K patt. that can be 


obtained, 





= 


VII, COMPARISON OF RESULTS 


Now to compare the results we have obtained, equation 5. 31, 
the spiral thrust time, is plotted against the energy added, equation 
9.22, infig. 3. Using a value of R=0.8 ,» which corresponds to 
a maximum ©,= 0.944 by equation 6.3, the restricteduperiace 
thrust time is computed using equations 6.2, 6.5, and 6.6 for values 
of eccentricity from zero to 0.95, with r = 0,95. Ehecorre- 
sponding increase of energy is computed using equation 6.7, and 
Bgese results are plotted in fig. 3. Since the total timiesmsust peseaqua 
for either method at a specified energy addition, we may compute 
equations 5.31 and 6.10 at specific energy levels and take their ratio 
to plot (€,/€,) in fig. 3. This is the ratio of acceleration for 
the discontinuous to the spiral that is necessary to accomplish the 
mission. 

It seems clear from fig. 3, in view of the acceleration ratios 
required for the mission and equations 2.13, 2.22, and 6, Saataat ie 
will be greater for the discontinuous method if ,K is the same for 
both methods and starting from the circular orbit with os = 0. 95% 
For these conditions, the discontinuous thrust method will provide 
less payload. Further, if we decrease R » the initial cutoff angle, 
we see from equation 6.10 that © must be increased for the same 
total time; thus, y* will be even Laucemror the ciscontinuous 
method, resulting in even less payload. 

Using equation 7.10 and fig. 3, we have evaluated equation 


7.15 for three diffexent final energies, expressed as eccentricity, for 
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various values of ©. These results are plotted in fig. 4, 

From reference (3), we find that a typical value of K would 
be about 150 pounds per kilowatt. This value will decrease as space 
power plant technology progresses. In reference (3), we also find 
K wees of 75 pounds per kilowatt-hour for nickel-cadmium batteries, 
which have very reliable charge and discharge characteristics, and 
14 pounds per kilowatt-hour for silver-zinc batteries, which are not 
as reliable for a large number of recharging cycles. From equation 


7,3, we see that both values of K must be divided by 


batt. 
/ 
3600 (4./R,) : to have the proper non-dimensional time units. 


Thus, we have the values for of 16.8 for nickel-cadmium 


K yatt. 
and 3. 13 for silver-zinc. 

We cannot improve over the spiral thrust programimewitn ciemes 
battery supply if the final energy must correspond to an eccentricity 
of 0.95, since in fig. 4 the maximum r Rye spall OLO1E 2 

Using equations 5.35 and 6.9, we compute the non-dimension- 
al thrust time for one cycle, corresponding to an initial cutoff angle 
eo: 0, 5 and e,=O 
Pre Ge dt 2 oe = ee (8. 1) 

Zl ee! (0.95) 
Using this and equation 7.14, we compute the values of Mi for the 
nickel-cadmium battery, for «= 7.5 and o&=10.0, respectively: 
p= l. to ) w= 1,93 (8. 2) 
and for the silver-zinc battery: 


hM= 0.47 ; tae (S5-3)) 
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From fig. 4, we read F Ry2=l. Zl] for the final eccentricity of 
eee and “X=7.5. Butirom equations (2.15550. 2Z,.200) ceo ee 
values of HH for either battery will require vy to be greater 
than Ys ; again, we find no advantage in the restricted perigee 


method of energy addition. However, for the final eccentricity of 0.7 
we read from fig. 4 rRye of 1.29 and 1939 for ex's) ch meno and 
10.0, respectively, From equations 7.15 and 8.3, we compute for 
this condition: 

pS oor Ss =} o10 (8. 4) 
d ng 


for the two values of co, respectively. From equation 2. 25, we can 
See that in this case the restricted perigee method has a slightly 
larger payload than the spiral method for this particular final energy. 
To see this result more clearly, we form the ratio of equation 2, 25 


for the two methods using 8. 4: 
(Mi+Mslg — (f— Yq)" 
UM te ioe (1- Loyy)* 


toe. = 10,0. Now, clearly, the numerator 0:0, 51s greater than 


(3,5) 


the denominator, so the ratio must be greater than one. 

Although we have found that the restricted perigee method of 
energy addition, using silver-zinc batteries, yields a slightly larger 
payload than the spiral method for certain values of final energy, we 
still have the problem of battery reliability. There is no question of 
the recharging reliability of the nickel-cadmium battery, but we 


found the specific mass, to be too high to be useful. the 


i Bocce 


silver-zinc battery could be used for a slight payload advantage, be- 





sae 
eause ol its lower specific mass, butits recharging teliaviiey ia 
imost be improved. 

We must finally compare the two methods of energy addition 
when reliable, rechargeable batteries are part of the payload but are 
available to be used in the thrust program. In this case, H=0 and 
we see immediately from fig. 4 that there are certain values of 
and final energy where the restricted perigee has a decided advantage 
eer the spiral method. 

As an example, let us determine the advantage where a 


nickel-cadmium battery, = 16.8, is part of the payload. 


Kpatt, 
Using equations 7.14 and 8.1, we compute values for (My ee /™4) 


for “™=7.5 and oc = 10.0, respectively: 


_M batt. eee M batt. _ 1.93 (8. 6) 


Ms Mg 
Even though the ratio of battery to power supply mass is rather 
Mage, it is not unrealistic since we could easily have a sateiiie 
where both battery and power supply are but a small fraction of the 
initial mass. From fig. 4 we read an y2= Foo fOr $e tase 


eccentricity of 0.4 and o = 7.5, and since =O 
% 
Ve 


Forming the ratio of equation 2.25 for the two methods, using 8.7: 


CMe Msi (=¥a)o (3. 8) 


(M,tMsle (1-14.10 ¥4)° 


= 1.10 (8. 7) 


Similarly, from fig. 4we read FR yz of 1.29 and 1. 39 for the 


final eccentricity of 0.7 and «'s of 7.5 and 10.0, respectively. 
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Since HM=0O  , we have immediately: 


Xs Xs 7 
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for the two values of « respectively. Again forming the ratio of 


equation 2. 25 for the two methods, using 8.9, we have 


(EE a eral (Mi+Ms)4 (- ¢4)° 


—w 
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for the two values of & =%7.5 and c= 10.0, respectively. The 
ratios in equations 8.8 and 8. 10 will depend on the value determined 
for ¥ , computed from equation 7.12, but they do show a definite 
advantage for the restricted perigee method of energy addition. 

For example, suppose the value for Yy is 0.25, that is, 
let 25 per cent of M, be propellant. Substituting in equation 8. 10, 
we compute ratios of 1.097 and 1. 131 for the two values of ™ , re- 
spectively. Now, this is about a 10 per cent gain in payload and 
structure mass by using the discontinuous rather than the continuous 
method of thrust programming, where the batteries are part of the 
payload. 

It seems rather surprising that we have found no advantage 
for the restricted perigee method of energy addition for the final ec- 
centricity of 0.95, where we found a definite advantage for the final 
eccentricity of 0.7. There are two reasons for this result. First, 
we can see from equation 6.9 that the travel time for 2/7 radians 
increases as the eccentricity increases. Thus, if we thrust for the 


same amount of time during each cycle, there 15 a greater percentage 





cA 
of coasting time as the eccentricity increases. This requires an in- 
crease of Cones | » as shown in tig. 3, te aceconiplisi ime 
mission, and in turn, we find a decrezse in the ratio of Galves, 
so that there is no advantage when @r reaches 0.95. Second, by 
restricting the cutoff angle to one radian or less, we do not take full 
advantage of the greater velocities provided by the elliptic orbit as 
fre eccentricity approaches one. this can be seen cleacl im tae 
where,at a cutoff angle of one radian and @ =0.95, the specific 
power of the restricted perigee is many times greater than the spiral 
fr1nod of energy addition, To obtain the tull advantage of this much 
Boearer Specific power, we must go to larger cutoit anplles: 

A study of the restricted perigee method of energy addition 
at larger cutoif angles would require the use of numerical methods. 
This would be rather involved, using equation 4.13 to evaluate AC 
for various values of ©, and € , then using this result fora 
corresponding evaluation of the change in energy, equation 4. 12, 

fter tabulating these ''cycle'' results, summations of time and energy 

could be made with a corresponding calculation of y* for various 
values of © . By making such a numerical study, we can clearly 
define the area where the restricted perigee method, with batteries, 
has an advantage over the spiral method of energy addition. We will 
also be able to determine if there is any area where an advantage ex- 


ists without the use of batteries. 





a ee 


Is. CONCEUaIoOn 


We have found that the discontinuous method of energy addi- 
ion, using a restricted perigee, has no payload ancesteGe mtr enmiacs 
ratio advantage over the spiral method when a specified energy is 
tiaparted in a specified time and the cutoii angle is eeatrie teas omen 
radian or less. There is the possibility of a very small advantage, 
fOr a certain range of final energies, if silver-2ingsoame tt een ee 
eluded in the power supply to store energy during the ceastmeponloac, 
moeoviding their reliability can be improved te amasaceepe iy se me a 
the other hand, a moderate decrease of the power-plant specific mass 
will rule out any advantage of including batteries in the power supply. 

In the case where batteries are included in the paylcadeesay, 
Bo periodically energize a transmitter after the satellite bas been 
placed in final orbit, we can definitely increase the payload and 
structure mass ratio for certain vaiues of final energy by using these 
batteries in the restricted perigee thrust program. tWhis does not ap- 
ply in general, obviously, since the mass fraction of the batteries re- 
quired in the payload must be compatible with the mass fraction of 
batteries required for the thrust program to accomplish the mission. 
When we do have a compatible battery requirement, then we can in- 
crease the payload and structure mass ratio by using the restricted 
perigee method of energy addition. 

Since we are restricting our investigation to low thrust de- 
vices, a numerical study of the restricted perigee method of energy 


addition at cutoff angles greater than one radian could be made using 





2Aae 

the perturbation equations developed earlier. Such a study would 
follow the procedure briefly outlined imethe last paracgrapheor cic pre 
seeding section, and is necessary to clearly define the area wcrc sine 


restricted perigee method can be used to advantage. 
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